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1.  Introduction. 


Let  X,  Y  be  Banach  spaces  and  A  C  R"  a  bounded  interval.  Let  F  ;  A  x  X  — *  K  be  a  smooth 
operate!.  The  nonlinear  equation 

(1.1)  F(A,u)  =  0, 

with  parameters  A  €  A  is  called  parametrized  nonlinear  equations. 

In  this  paper  we  deal  with  the  parametrized  nonlinear  equation  F  :  A  x  Hq{J)  — ► 
with  one  parameter  A  €  A  C  R  defined  by 

(1.2)  F(A,u)  =  0,  (A,u)  €  A  X  Fo^J), 

(1.3)  <  F(A,u),v  >:=  j{a(X,x,u'(x))v'(x)  +  f(X,x,u{x))vix)]dx,  Yv  €  HliJ), 

where  J  :=  (6,  c)  C  R  is  a  bounded  interval,  a, /:AxJxR— 'R  are  sufficiently  smooth 
functions,  and  <  •,  •  >  is  the  duality  pair  of  and  Hq(J).  Since  F  is  a  second  order 

differential  operator  in  divergence  form,  finite  element  solutions  of  (1.2)  are  defined  in  a  natural 
way. 

In  [TBl],  the  above  problem  was  concerned,  and  a  priori  error  estimates  of  finite  element 
solutions  were  established.  In  this  papar  we  try  to  develop  a  postriori  error  estimates  of  the 
finite  element  solutions  of  (1.2)  and  (1.3)  using  a  priori  estimates  obtained  in  [TBl]. 

The  basic  idea  is  as  follows.  Suppose  that  we  want  to  solve  the  nonlinear  equation 

(1.4)  <K{x),v>=0  for  'dv  e  Hq(J), 

O 

where  K  :  Hq{J)  — ►  is  a  smooth  nonlinear  operator.  Let  S\  C  Hl{J)  be  a  finite  element 

P 

space  and  €  5/,  the  finite  element  solution,  that  is, 

(1.5)  <  K{xh),Vh  >=  0  for  6  5/,. 

Then,  we  consider  the  linearized  equation 

(1.6)  <  DK{xh)^,v  >=  -  <  K{xh),v  >  'iv  e  Hq(J), 

where  DK(xh)  is  the  Frechet  derivative  at  which  is  assumed  to  be  an  isomorphism  between 
H^{J)  and  R-HJ). 

Let  X  6  Hq(J)  be  the  exact  solution  of  (1.4).  In  Section  3  we  will  see  that  the  magnitude 
(IV'II  represents  the  error  (|i  -  r^(|,  that  is, 

Ik  -  2:a||  <  \\rp\\{l  +  o(l)). 
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Of  course,  the  exact  solution  i/  6  Hq{J)  of  (1.6)  should  be  approximated  by  a  certain  way 
in  general.  We  will  consider  a  finite  element  solution  of  (1.6).  We  observe,  however,  that  the 

O 

finite  element  solution  rph  of  (1.6)  over  Sk  defined  by 

O 

<  DK{xh)‘iph,^h  >=  -  <  K{xh),VH  >,  Vuh  6  Sk 

is  a  zero  function  because  of  (1.5).  Therefore,  estimating  the  magnitude  ||i/'||  is  equivalent  to 
estimating  ||0  —  the  error  of  the  finite  element  solution  i>k- 

Hence,  if  we  have  certain  methodology  for  a  posteriori  error  estimates  for  the  linearized 
equation  (1.6),  we  have  a  posteriori  error  estimates  for  the  original  nonlinear  equation  (1.4).  In 
a  short  sentence,  the  principle  obtained  here  is  that 

"If  we  have  a  posteriori  estimates  of  linear  equations, 
we  have  a  posteriori  estimates  of  nonlinear  equations.” 

Let  (A,u)  be  the  exact  solution  of  (1.2)  and  (A/,,Uh)  be  the  finite  element  solution  corre¬ 
sponding  to  (A,u).  Usually,  it  is  observed  that  the  error  |A  -  A^j  is  much  smaller  than  the  error 
||u  —  Ufcll.  In  Section  4  we  obtain  elaborate  error  estimates  of  |A  -  A/,|  which  verify  the  above 
observation. 

In  Section  5  practical  aspects  of  our  a  posteriori  estimates  and  some  numerical  examples 
are  given.  In  the  computation  of  our  numerical  examples  the  continuation  program  package 
PITCON  (see  [R])  developed  by  Rheinboldt  and  his  colleagues  is  used. 

This  paper  is  a  revision  of  a  part  of  one  of  the  authors  Ph.D.  dissertation  [T]. 

2.  Assiunptions  and  A  Priori  Estimates. 

In  this  section  we  sumerize  the  results  obtained  in  [TBl].  Throughout  this  paper,  we  use  same 
notation  as  in  [TBl]. 

Here,  we  deal  with  the  nonlinear  operator  F  ;  A  x  Wq’°°  —  *  by,  for  A  6  \  and 

(2.1)  <  F(A,u),u  >:=  j^[a{X,x,u\x))v\x)-k-  f{\,x,u{i))v(x)\dx,  'iv  €  ^ 

where  <  •,  •  >  is  the  duality  pairing  between  and  Wg  *.  Then,  our  problem  is 
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Problem  2.1.  Solve  the  following  equation:  Find  A  €  A  and  u  6  such  that 

<  r(\,u),v  >=0,  VvelVo^’^.  □ 

For  F  being  well-defined  and  smooth  we  require  several  conditions  to  a  and  /.  Let  a  = 
(011,012)  be  usual  multiple  index  with  respect  to  A  and  y.  That  is,  for  or  =  (01,02),  D°a(\,x,y) 
means  Qy°‘t  ^(^,  v)- 

Let  d  >  1  be  an  integer.  For  o,  |o(  <  d,  we  define  the  maps  A“(u)  and  F“(u)  for  u  € 


by 

(2.2) 

A"(u)(i) 

:=  I>“a(A,x,u'(x)), 

(2.3) 

F“(u)(x) 

:=  D“/(A,x,u(x)). 

We  then  assume  that 

Assumption  2.2.  Foi  all  o,  |o|  <  d,  we  suppose  that 

(1)  For  almost  all  x  e  J,  D‘‘a(X,x,y)  and  D°‘ f{\,x,y)  exist  at  all  (A,y)  6  A  x  R.  and  they  are 
Caratbeodory  continuous. 

(2)  The  mapping  A“  defined  by  (2.2)  is  a  continuous  operator  from  Wq’°°  to  and  the  image 
f\‘^(U)  C  L°°  of  any  bounded  subset  U  Q  \  x  Wq’°°  is  bounded. 

(3)  The  mapping  F"  defined  by  (2.3)  is  a  continuous  operator  from  to  ,  and  the  image 

F°(U)  C  of  any  bounded  subset  U  C  \  x  is  bounded.  □ 

We  define  the  subset  5  C  A  x  Wo'*”  by 
(2.4)  5  :=  {(A,u)  €  A  x  a,(A,  x,  u'(x))-»  €  I*}. 

Since  the  mapping  A  x  Wq'°°  3  (A,u)  1—  aj(A, x, u'(x))  6  L°°  is  continuous,  we  have 

Lemma  2.3.  If  a  and  f  satisfy  Assumption  2.2  with  d>  1,  S  is  an  open  set  in  A  x  U’o 

□ 


From  the  standard  theory  of  Fredholm  operators,  we  obtain  the  following  theorem; 

Theorem  2.4.  Suppose  that  a  and  f  satisfy  Assumption  2.2  with  d  >  1.  Then  in  S.  the 
operator  F  :  5  — •  defined  by  (2.1)  is  a  nonlinear  Fredholm  operator  of  index  1.  □ 
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We  define  the  subset  7l(F,S)  C  S  hy 


(2.5) 


Tl(F,S)  :=  {(A,u)  €  5|Z)i’(A,u)  is  onto}. 


The  elements  of  5)  and  f('/J(J’,5))  are  called  regular  points  and  regular  values, 
respectively.  By  Theorem  2.4,  we  can  apply  the  Fink-Rheinboldt  theory  ([FRl],[FR2],[R])  to 
the  operator  F  and  obtain  the  following. 

Theorem  2.5.  Suppose  thit  a  a.nd  f  satisfy  Assumption  2.2  with  d  >  1.  Let  e  €  FiTli  F,S)). 
Then 

M  =  M^  :=  {(A,u)  e'fl(F,S)  |  F(A,u)  =  e} 

is  a  oae-dimeasionaJ  C^-manifoId  without  boundary.  Moreover,  for  each  (A,  u)  6  M.  the  tangent 
space  T(x,u)A4  at  (A,u)  is  KetDF{X,u). 

Therefore,  if  0  €  F{1l(F,S)),  the  solutions  of  Problem  2.1  form  a  one-dimensional  C*^- 
manifold  without  boundary  in  7l(F,S).  □ 

In  the  sequel  of  this  paper  we  always  assume  that  0  6  F(1Z(F,S)),  that  is,  .Vto  ^  0. 

For  the  regularigy  of  (A,u)  6  Mq,  we  need  additional  assumptions.  Let  p’,  2  <  p"  <  x  be 
taken  and  fixed. 

Assumption  2.6.  Under  Assumption  2.2  with  d  >  1,  we  assume  that 

(1)  For  aii  A  €  A,  the  functions  a(A,  •,  •),  Oy(A,  •,  •) :  7  x  R  —  R  are  continuous. 

(2)  For  all  (A,y)  e  A  x  R,  there  exist  ai(A,i,y)  for  almost  all  x  £  J  and  are  Caratheodory 
continuous. 

(3)  The  composition  functions  /(X,x,u(x)),  ax(X,x,u'(x))  are  in  L^’  for  any(X,u)  €  A  x  Wq 
Moreover,  for  any  bounded  subsets  K  C  A  x  Wg  '*, 

{f(X,x,u{x))£L^’\(X,u)eK},  {MA,x,a'(z))e  i'’’!  (A,u)€  A'} 

are  bounded  in  L^' .  □ 

Theorem  2.7.  Under  Assumption  2.2  and  2.6,  we  have  u  G  1^2, p*  (A.u)  G  .Vto- 

Moreover,  for  all  bounded  closed  subsets  M  C  Mq,  there  exists  a  constant  A'(.U)  such  that 
sup_||u||wri,p*  <  A'(M). 
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Let  Sh  C  Hq  he  a.  finite  element  space.  We  define  the  finite  element  solutions  of  Problem  2.1 
by 

.  • 

Problem  2.8.  Find  A;,  G  A  and  6  5/,  sack  that 

<  >=  0,  'ivheSn-  □ 

Let  cr  €  L°°  be  such  that  a(i)  >  c  >  0  for  all  i  6  J  Let  (•,  la  be  the  inner  product  of  Hq 
defined  by  {u,v)a  '■=  j  au'v'dx  for  u,v  €  Hq.  Define  the  isomorphism  Ta  € 
by  <  Ti,v  >=  {TQV,'<j)a,  Vv  G  for  rj  G  Also,  define  the  canonical  projection 

-♦  Sh  by  {ip  -  =  0,  €  Sh  for  ip  G  H^.  Then,  we  observe  that,  for  any 

O 

Vh  €  Sh  and  any  v  ^  Hq, 

<  F{\h,Uh),Vh  >=  0  <=><  T~^Tl'jlTaF(Xh,yh),v  >=  0- 
Following  the  Fink-Rheinboldt  theory  we  define  7^  :  A  x  by 

ThiK u)  :=  (/  -  (A,  u), 

where  I  is  the  identity  of  and  :=  T“^II^To. 

Lemma  2.9  ([R, Lemma  5.1]).  The  operator  F°  satisfies  the  following: 

(1)  f^(A, «)  =  0  for  some  (A,  u)  G  A  x  Pg  i/  and  only  if  (A,  u)  G  A  x  Sh  and  F;,(  A.  u)  =  0. 

(2)  7^  is  a  Fredholm  operator  of  index  1  on  \  x  Hq.  □ 

By  Lemma  2.9,  we  have  the  following  theorem  as  a  consequence  of  the  Fink-Rheinboldt 
theory. 

Theorem  2.10.  Suppose  that  F  is  C**  mapping  (d  >  1).  Then  the  set  of  the  finite  elements 
solutions  of  Problem  2.8, 

Mh  :=  {(Afc,uO  G  Tl{Fh,\  x  H^)) \Fh{Xh,Uh)  =  o}  . 
is  aC^  manifold  without  boundary.  □ 

For  a  priori  error  estimates  of  the  finite  element  solutions,  we  always  assume  the  following. 
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Assumption  2.11.  We  assume  tha.t 

(1)  Assamptioa  2.2  with  d  (i.e.  F  is  Fredholm  map). 

(2) 0e  F(7l(F,S))  (i.e.  Mo  9^  ^). 

(3)  Assumption  2.6  (i.e.  u  6  W^'^' ,  2  <  p’  <  00  for  any  (A,  u)  6  Mq). 

(4)  Sh  is  tegular  and  Um  inf  llu  —  Uhllui  =  0,  for  any  u  €  Hi. 

o 

(5)  The  triangulation  of  5/,  (in  one  dimensional  case,  the  partition  of  J  into  small  intervals) 
satisfies  the  inverse  assumption  [C,pl40].  □ 

In  the  sequel,  we  denote  by  ft/,  :  —  S/,  the  interpoiant  projection. 

Theorem  2.12.  Suppose  that  Assumption  2.11  holds  for  d  >  2.  Also,  suppose  that  Mq  C 
Mq  is  a  compact  regular  branch,  that  is,  there  is  a  compact  interval  A  C  A  and  map 
A  9  A  t-»  u(A)  6  Wq'°°  such  that 

Mq  —  |(A,u(A))  6  Mq  I  DuF(\,u(\))  is  an  isomorphism  for  VA  €  a|  . 

Then,  there  exists  the  corresponding  finite  element  solution  branch  Mh  C  Mh  which  is 
parametrized  by  the  same  A  €  A  and 

||ft/.u(A)-«„(A)|ljy.  <  Kohh\ 

IKA)-«/.(A)||^.  <  A-i|MA)-ft„u(A)||^., 

MX)  -  u^{X)\\^i.oo  <  K2h^ 

for  all  X  €  A,  w(A)  €  Mq,  u/,(A)  €  Mh,  and  77  with  0  <  r/  <  Here,  Kq,K\,K2  >  0  are 
constants  independent  of  h  and  X. 

Moreover,  we  have 

MhCn(F,S).  □ 

Theorem  2.13.  Suppose  that  Assumption  2.11  holds  for  d  >  2.  Let  Mq  C  Mq  be  a 
connected  compact  subset  with  the  following  properties: 

(1)  DxF(X,u)  ^  0  for  any  (X,u)  €  A<o. 

(2)  There  exist  xq  €  J  such  that  DG(X,u)  defined  by,  for  given  7  €  R, 

(2.6)  G(X,u)  :=(u(xo)-y,F(X,u)),  (A,  u)  6  5 

(2.7)  DG(X,u)(t,^i>)  =  (^i>(xQ),DF(X,u)(t,ii>)),  t  €  R.  €  H’q 
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is  aa  isomorphism  at  aii  (A,w)  €  A<o- 

Then  A(o  Js  pa.ra.mettized  by  y  =  u(xo).  We  assume  without  loss  of  generality  that  the 

Q 

above  xq  is  a  nodal  point  of  Sh  for  all  sufficiently  small  h  >  0. 

Then  there  exists  the  corresponding  Unite  element  solution  branch  Mf,  C  Mh  which  is 
parametrized  by  the  same  y,  that  is,  U|,(7)(io)  =  7  and  T),(Afc(7),  u/,(7))  =  0  for  any  y. 
Moreover,  we  have 

|A(7)  -  Afc(7)|  +  ||n;.«(7)  -  «k(7)||//.  < 

|A(7)  -  Afc(7)|  +  ||u(7)  -  «fc(7)||//^  <  A'4||w(7)  -  n/,u(7)|lH„>  - 
|A(7)  -  A/i(7)|  +  ||u(7)  -  w/,(7)||vv,j.oo  <  A'sh’', 

Mh  C  nF^S), 

for  all  y  =  u(xo),  (X(y),u(y))  €  A<o.(A/,(7).«a(7))  €  Mt,,  and  r?  with  0  <  tj  <  \.  Here. 
A'a,  K^,  A's  are  positive  constants  independent  of  h  and  y.  □ 

3.  A  Posteriori  Error  Estimates. 

In  this  section  we  consider  a  posteriori  error  estimates.  Before  going  into  our  problem,  we 
observe  an  error  estimate  in  a  general  Banach  space  setting. 

Let  X  and  Y  be  Banach  spaces  and  V  Q  X  open.  We  consider  a  generic  C-  mapping 
K  :  V  -*  Y  such  that  DK{x)  €  C{X,Y)  is  an  isomorphism  at  each  x  €  V^,  and  D~I\(x)  is 
bounded  on  bounded  subsets  in  V . 

Suppose  that  we  are  considering  the  equation 

(3.1)  K{t)  =  0,  t^V. 

Let  tgx  €  V'  be  an  exact  solution,  i.e.  A"(t£x)  =  0,  and  t^p  €  V'  an  approximate  solution,  i.e. 
K{tAp)  «  0.  Note  that,  since  DK{tEx)  is  an  isomorphism,  tgx  €  V'  is  isolated,  that  is.  there 
is  no  other  solution  of  (3.1)  in  the  small  enough  neighborhood  of  tgx- 
From  elementary  calculus  on  Banach  spaces,  we  have 

(3.2)  0  =  A  (t^p)  +  Dh{tAp)z  +  -  (1  -  s)‘ K{tAP  +  sz)ds^  (;.::). 

where  z  tgx  —  t^p-  Let  us  consider  the  following  linearized  equation: 

(3.3)  0=  K{tAP)  +  DK{tAP)z,  5€l/. 
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From  (3.2)  and  (3.3)  we  obtain 


(3.4)  DK{tAP){z  “  ~  ~  +  32)^5^  (z,  2), 

(3.5)  2  -  2  = -^I>A'(f^/>)~^  (1  -  s)^D^ii'(t^/>  +  3z)ds^  (2,2)  , 

and 

11^  -  51U  <  ^IPA'(l^PrMl£(K.X)  [  \\D^K{tAP  +  32)ll£(xxx,y)d5ll2|l2;,. 

Since  D^K  €  C(X  x  X,  Y)  is  bounded  on  bounded  subsets,  there  is  a  constant  M  such  that 
II2  —  2|lx  <  ^lUllxi  3.®^  obtain 

(3.6)  ll^llx<||i||x(l  +  0(||2||x)). 

By  the  argument  in  [TBl, Section  4],  we  know  that,  at  each  (A,u)  6  Mq,  we  have  either 
Case  1:  Keri)uF'(A,  u)  =  {0}  and  D\F{\,u)  €  ImZ)uF(A,u),  or 
Case  2;  dimKerJPuf’(A,  u)  =  1  and  D\F{\,u)  ^  ImD„i'(A,u). 

Now,  let  us  suppose  that  we  are  in  Case  1.  To  apply  (3.6)  we  set  up  the  following 
X  =  Y  = 

X  dV  =  niF,S)  defined  by  (2.5), 

K{u)  =  F(A,ti)  for  given  and  fixed  A  €  A. 

From  (3.6),  we  have 

!|U  —  U;,||^,J,oo  <  ||f^|ly(/^l.oo(l  +  C(||f71|j^j^J.oo)), 
where  (A,ti)  €  A4oi  (A,ti/,)  €  and  U  is  the  exact  solution  of  the  linearized  equation 

(3.7)  oo<  .C>„F’(A,u^)f7,u  >1= -«< /’(A,Uh),  V  >1,  Vv  6  I'V’g  ^ 

By  Theorem  2.12,  we  have  lim  F(X,Uh)  =  0  and  lim  ||f7||,yi.oo  =  0.  Therefore,  we  obtain 

Theorem  3.1.  Suppose  that  (A,u)  6  A4o  is  on  a  regular  branch.  Then  we  have,  under 
Assumption  2.11  with  d  >2, 

||u  —  U^II^J.oo  <  ||(7|l(^^l,oo(l  +  0(1)).  □ 

We  next  consider  an  a  posteriori  error  estimate  in  /Tg-norm.  From  (3.5).  we  have 

2  _  =  -iD„F(A,  u,)-'  [(^\l  -  5)-D;„F’(A,  +  sz)ds^  (2.  2)  . 
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where  z  :=  u  -  Uh-  Recall  that  £>u/’(A,u)  6  C{Hq,  H~^)  is  an  isomorphism,  and 

oo<  (1  -  sfDl^F{\,Uh  +  sz)d^  {z,z),v  >1 

=  j {aXx)z^v'  +  0i{x)z^v)dx,  'iv  € 
fX  y  1 

where  ai(r)  :=  /  (1  -  s)^ayy{X,x,u’f^  +  sz')ds  and  0\{x)  :=  /  (1  -  s)^/yy(A,  x,  u/,  +  sz)ds. 
Jo  Jo 

Hence,  we  get 

(^l\l-sfDl,F(X,u^,  +  sz)dsyz,z)  <  ||ai|U=o||z'2|L. +  ||AIL.||-ili» 

and 

11^  ~  ^  C’2ll~llw^'=o|kll^ij  • 

Therefore,  by  Theorem  (2.12),  we  obtain  ||  <  ||C/^||^^  (l+Ca/i'')  for  any  r]  with  0  <  r?  <  4. 

and 

Theorem  3.2.  Suppose  tha,t  (A,u)  €  ^«o  is  on  a  regular  branch.  Then  we  have,  under 
Assumption  2.11  with  d>2, 

11^  ~  “hlltfij  <  +  o(l))-  O 

Next,  let  (A,u)  €  AAq  such  that  D\F{\,vl)  ^  0.  Then,  by  Theorem  2.13,  there  exists  a 

O 

nodal  point  xq  ^  J  o[  such  that  the  Frechet  derivative  DG(\,u)  of  the  mapping  G’(  A.  u)  := 
(u(xo)  -  7.  F{X,u))  is  an  isomorphism  of  R  x  to  R  x 

We  consider  the  following  problem. 

Problem  3.3,  For  given  y  €  R,  find  u  €  Wj  *  and  A  €  A  such  that 

oo<  •f’(A,  u),  t;  >1=  0,  Vu  6  Wo’\  and  u(xo)  =  7.  □ 

Problem  3.3  'orresponds  to  the  equation  G(A,u)  =  (0,0).  Naturally,  we  define  the  finite 
element  solution  for  Problem  3.3  by 

0 

Problem  3.3f  £.  For  given  7  €  R,  find  u/,  e  S^  And  Xh  6  A  such  that 

O 

<  F(Xh,Uk),Vh  >=  0,  Vu/,  €  5(,,  and  Uh(xo)  =  7.  □ 
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Since  DG{\,u)  is  an  isomorphism,  we  can  apply  (3.6)  to  Problem  3.3  and  3.3yr£.  The 
linearized  equation  is 

(3.8)  {U{xq),9DxF{\h,uh)+  D^F{XH,^h)U)^iO,-F{XH,Uh)),  9eR,  UeW^'^. 

It  follows  from  Theorem  2.13  that  lim  f’(A/,,u/,)  =  0  and  lim(||17||^i,<x)  +  |0|)  =  0.  We  set 

A— »0  A-»0  '0 

X  =  Rx  Y  =  Rx 

K{X,u)  =  G{X,  u)  for  given  7  €  R 
By  (3.6),  we  have 

Theorem  3.4.  Suppose  tha.t  (A,u)  G  Mo  satisfies  D\F{X,u)  ^  0.  Then  we  have,  under 
Assumption  2.11  with  d  >2, 

lA  -  Afcl  +  lln  -  <  (\9\  +  (1  +  0(1)).  □ 

We  can  get  an  a  posteriori  error  estimate  in  Jg-norm.  Rewriting  (3.4)  in  the  above  setting, 
we  have 

(3.9)  DGiXk,UH)it  “  J  “  sfD^GiXh  +  st,UH  +  3z)ds^  («.;)'. 

where  t  :=  A  -  A/,,  2  :=  u  —  u/,,  and 

(3.10)  D^G{X,u)  =  i0,D^F{X,u)). 

Since 

<  D'^F{X,u){t,  z)^,v  >  =  jjaxx{X,x,u')v'  +  fxx{X,x,u)v]dx 

+  2t  j^[axy{X,  X,  u')z'v'  +  fxy{X,  x,  u)zv]dx 
+  J^[ayy(X,  X,  u')z'^v'  +  fyy{X,  X,  u)z-v]dx, 

we  easily  obtain 

(3.11)  ||Z)2F(A,,  +  St,Uh  +  52)(t,2)2l|//-l  < /l|«p  +  B|t|l|2||^^  +  Ci|2|l^J.=o||-'||/^l 

for  any  s  6  [0, 1],  where  A,  B,  C  are  constants  independent  of  h. 

Therefore,  from  (3.9)-(3.11),  there  exists  a  constant  M  such  that 

It  _  0|  +  11^  _  <  A/(|^|2  +  (1^1  +  ||,||^.,^)||.||^,  ), 

and  we  obtain 
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Theorem  3.5.  Suppose  tha.t  (A,u)  €  Mq  satisfies  DxF{X,u)  ^  0.  Then  we  have,  under 
Assumption  2.11  with  d  >2, 

-  '“hlUo*  -  (l^l  +  II^IIko')  (1  +  ° 

Now,  to  get  another  a  posteriori  estimate,  we  consider  the  following  auxiliary  equation:  find 
W,f  €  Hq  such  that 

^  ^  r  -iaHix)W:,y  +  =  -Fh  -  nKn  on  J  -  {xo}, 

(3.12)  ( 

I  lV,(xo)  =  0, 

where  0^(1)  :=  a5(Afc,  x,  ti;,(x)),  /3/,(x)  :=  /y(Afc,  i,  Ufc(x)),  and 

Fh  :=  -a(A;„x,tifc(x))' + /(A;,,x,u^(i)), 

Kh  ■=  -ai(Ah,i,uUx))' + /A(Ah,x,u;,(x)). 

Since  (3.12)  is  equivalent  to  DG{\h,Uh){Q,  IV,)  =  (0,  —Fh—rjKhO)  we  see  that,  for  sufficiently 
small  h  >  0,  (3.12)  has  a  unique  solution  W,,  for  7?  6  R  which  is  sufficiently  close  to  9  (see  the 
proof  of  Lemma  4.3). 

Note  that,  even  if  u'f^{x)  is  not  continuous, 

$,,(x)  :=  afc(x)iy,J(x)  +  77aA(Aft,x,u{,(i))  +  a(A/,,x,u),(i)) 

is  continuous  on  7  —  {iq}-  Then,  we  define  the  ‘jump’  J{t})  at  x  =  xq  by 

(3.13)  7(77):=  lim  #,(i)-  lim  ^Jx). 

X'^xo  + 

From  (3.8)  and  (3.12),  we  clearly  have  7(0)  =  0.  Let  U  :=  Wq.  Then,  we  claim  that 
Lemma  3.6.  We  have 

1^1  +  ^  ll^llvv'^‘«>(l  +  o(l))i 

M  +  \\U\\Hi  <  l|f^||^.(l+0(l)). 

Proof.  Integrating  (3.12)  by  part,  we  have 

(3.14)  <  DuF{Xh,'^h)U,v  >=  -  <  F(Aft,u/,),u  >  +i7(xo)7(0), 

O 

for  any  v  ^  Hq.  Taking  €  5/,  with  '>p{xo)  =  1  and  fixing  it,  we  obtain 

(3.15)  <  D„F(\H,Uh)U,i^  >=  J{0). 
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Let  2  €  /To  be  the  function  which  satisfies  z{xq)  =  0  and 

<  DuF{Xh,UH)z,v  >=<  DuF(\h,Uh)rl^,v  >,  Vi;  €  ^obo], 
where  /fo[io]  '•=  €  /Tq  |  w(xo)  =  0}.  Since  U{xo)  =  0,  we  have 

(3.16)  <  DuF{Xh,^th)z,U  >=<  DuF{Xh,Uh)ip,U  >  . 

It  follows  from  (3.14)  that 

(3.17)  <  D^F{Xh,Uh)U,vuh>=  -  <  F{Xh,ut,),Wh>=0,  Vw;,  G  5/,[xo], 

where  5/,[io]  :=  jwfc  €  \  Vfc(io)  =  o|. 

Since  DnF{Xh,Uk)  is  self-adjoint,  we  obtain 

(3.18)  <  Z?„f’(A;„  Uh)(z  -  Wh),  U  >=<  DuF{Xh,Uf,)'ip,  U  > 
by  subtracting  (3.17)  from  (3.16).  Combining  (3.15)  and  (3.18),  we  obtain 

|J(0)|  =  I  <  DuF{XH,UH)iz  ~wh),U  >  I,  6  5fc[io]- 
Now,  letting  wh  be  the  finite  element  solution  of  z,  we  get 

k(0)l  <  C'lll^  -  bm  Ik  -  =  0, 

where  Ci  is  a  constant  independent  of  h. 

Since  (3.12)  is  a  linear  equation  with  respect  to  Wr,,  the  implicit  mapping  rj  Wj,  defined 
by  (3.12)  is  C°°  and  therefore  there  exists  the  derivative  of  W,^  with  respect  to  tj.  We  denote  it 
by  djfWr,.  We  see  that  dr,^,,  6  Hq,  and  it  satisfies 

(3  19)  !  =  -A\  on  J  -  {zo}, 

1  9,W,(io)  =  0, 

Note  that  integrating  (3.19)  by  part  shows  us  that  J(j7)  is  differentiable  and  we  have 
<  DxF{X^^,u^^)  +  DuFiXh,Uk)d,,W,„v  >=  v{xo)J'{t]),  'iv  G  Hq. 

We  remark  that  J'{rt)  yt  0.  If  J'{r])  =  0,  we  would  have 

{dr,Wr,{xo),  DxF{X,^,uh)  +  DuF{Xh,Uh)d,,W^)  =  (0,0). 

Since  DG{Xh,Uh)  is  an  isomorphism,  we  obtain  the  contradiction  {\,d,,W,j)  =  (0,0). 
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Therefore,  we  obtain  —^(0)  =  J{9)  —  J(0)  =  9J'{^0),  0  <  ^  <  1,  and 


(3.20)  1^1  <  \fmrvm  <  ^211.'  - 

On  the  other  hand,  we  have 

(3.21)  \\U  -  U\\y^rl«‘  <  ^3|^l, 

because  U  =  Wg,  U  =  Wq,  and  the  mapping  t)  •-*  W,,  is  C°°.  Combining  (3.20)  and  (3.21),  we 
conclude  that 

<  ||cri^.,cc(i  +  o(i)). 

The  second  inequality  is  proved  by  the  same  manner.  □ 


Theorem  3.7.  Suppose  that  (A,u)  6  Mq  satisfies  DxF{X,u)  ^  0.  Then  we  have,  under 
Assumption  2.11  with  d  >  2, 


|A  -  A/i|  4- ||u  -  <  ||t^||n,.j.oo(l  +  o(l)), 

+  11“  “  S  II^IUl(l  +  o(l)). 


where  17  €  Hq  is  the  exact  solution  of  the  following  equation: 

J^[ay{\h,x,u'f^)U'v'  +  fy{XH,x,Uk)Uv]dx  =  -  <  F{Xh,Uh),v  >,  Vt;  e  H^[xo], 
U{xo)  =  0.  □ 


(3.22) 


Remark  3.8.  From  Theorem  3.1,  3.2,  3.4,  3.5,  and  3.7,  a  posteriori  estimates  of  the  error 
||u  —  u^ill  and  |A  —  A/,|  +  ||«  —  u/,||  are  reduced  to  estimates  of  ||f^|l,  |0|,  and  ||^’||.  respectively. 
We  will  compute  the  finite  element  solutions  17*,  0*,  17*  of  the  equations  (3.7),  (3.8),  and  (3.22) 
and  see  ||17*||,  |5*|,  ||17*||  instead  of  l|t7||,  \e\,  ||t7||. 

We  must,  however,  notice  that  the  right-hand  sides  of  those  linearized  equations  (3.7). 
(3.8),  and  (3.22)  are  the  terms  F(A,u*)  and  F(A*,u*).  By  the  definition,  those  terms  vanish 

O  O 

for  u*  e  S*,  that  is,  <  F(A,u*),u*  >  =  <  F(A*,u*),t;*  >=  0  for  any  i;*  €  S*.  Hence,  the  finite 

O 

element  solutions  of  (3.7),  (3.8),  and  (3.22)  over  S*  would  be  just  the  zero  functions  and  they 
would  be  useless  to  estimate  ||17||,  \d\,  and  ||^||. 
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We  use  the  different  finite  element  space  5*  such  that  5/,  C  5*  to  avoid  this  difficulty.  That 
is,  to  compute  the  finite  element  solutions  of  (3.7),  (3.8),  and  (3.22),  we  use  the  refined  mesh 
or  higher  order  polynomials  on  each  finite  element.  Then,  we  will  obtain  nonzero  finite  element 
solutions  and  \\Uhi\\,  \\Uhi\\  indicate  the  error  |lu-u;,||  and  |A- Ah|  +  ||u-u/.||, 

respectively. 

Note  that,  in  the  computation  of  the  ((f4ill,  |^k,  |.  ||^h,  ||)  we  do  not  need  to  solve  the  entire 
problem.  The  estimations  of  ||£^a,  1|,  1,  ||^fc,ll  are  done  by  an  element- by-element  approach 

(see  [BR]). 

The  details  of  the  practical  computation  will  be  presented  in  Section  5.  □ 


4.  Elaborate  Error  Estimates  of  |A  —  A/,|. 

Sometimes,  one  may  want  to  estimate  only  the  error  |A-  A/,|.  Usually,  it  is  observed  that  |A  -  A^| 
is  much  smaller  than  ||u  —  u/i||^i.  In  this  section  we  develop  two  elaborate  error  estimates  of 
|A-Ak|. 

Let  (A,  ti)  6  Mq  be  such  that  (A,  a)  is  around  a  turning  point  or  on  a  ‘steep  slope',  that  is, 
DxF(X,  a)  ^  0.  Let  the  nodal  point  iq  €  J  of  5^  be  taken  so  that  DG{X,  a)  6  £(R  x  Wq R  x 
W^’*”)  is  an  isomorphism  (see  Theorem  2.13).  Let  7  :=  a(xo).  Then,  (A,  a)  is  a  solution  of  the 
following  problem: 


(4.1) 


Problem  4.1.  Find  a  €  **  A  €  A  such  that 

I  ^[a(A,i,u'(i))v'(z)-|-/(A,i,a(i))a(x)]di  =  0,  Va  6  \ 

u(io)  =  7. 


□ 


By  Theorem  2.13,  it  is  guaranteed  that,  for  sufficiently  small  h  >  0,  there  exists  a  locally 
unique  solution  (Ah,Ufc)  €  Mh  of  the  following  problem  around  (A,  a)  6  Mq. 


(4.2) 


Problem  4.1f£.  Find  a  €  5/,  and  A^  €  A  such  that 

J^[a(Xk,x,u'^(x))vx(x)  +  f(Xh,x,Uh(x})vh{x)]dx  =  0,  Va/,  6  5/,. 


«fc(io)  =  7-  □ 

To  estimate  the  error  |A  —  A/i|  we  introduce  the  following  auxiliary  equation. 
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Problem  4.2.  Find  u  €  such  that 

f  -(a(A/i,x,u'(x))y  + /(^^,x,u(x))  =  0  on  J  -  {xq}, 

(4.3)  { 

I  u(xo)  =  7-  ° 

On  the  existence  of  the  solution  u  of  Problem  4.2,  we  prove  the  following  lemma. 

Lemma  4.3.  Suppose  that  Assumption  2.11  holds  for  d  >  2.  Let  (Aq.uo)  €  Mo  be  such 

O 

that  DxF{Xo,uq)  ^  0.  By  Thorem  2.13  we  can  take  a  nodal  point  xq  e  J  of  S/,  such  that 
i?G( Ao,«o)  defined  by  (2.7)  is  an  isomorphism. 

Then,  there  exist  e  >  0  and  a  unique  map  (Aq  -  e,  Aq  +  e)  B  Ah—  u;(A)  €  *  which 

satisfies 

(4.4)  -  (a(A,i,u;(A)'(i)))'  + /(A,i,u/(A)(i))  =  0,  w(X)(xq)  =  y, 
where  y  ;=  uo(xo)- 

Prooj.  Let  W^o’^[io]  :=  G  *^o’^l^(2^o)  =  0}-  First,  we  note  that  u;(A)  satisfies  (4.4)  if  and 
only  if  it  is  the  solution  of  the  following  equation: 

J  ^[a(A,i,u/(A)'(i))u'(i)  + /(A,i,i£/(A)(i))v(x)](ii  =  0,  Vi;  € 

[  u/(A)(io)  =  y. 

Let  (W''o’^[io])’  be  the  dual  space  of  VFo’^[xo].  We  define  the  mapping  G  :  X  IFJ  *  — 
R  X  (Wq’^[io])’  by  G(X,w)  :=  ^u;(xo)  -  7, F’(A,u;)j^,  J.  Then  we  have,  for  u>  €  Wj  ”, 

(4.5)  Du,G{X,w)i>  =  ^^(xo),(DuF'(A,w)t(')||^,.,j^^J  • 

If  we  show  that  D^G{Xo,uo)  €  £(Wg’'”,R  x  (Wj’*[xo])*)  is  an  isomorphism.  Lemma  4.3  is 
proved  by  the  implicit  function  theorem. 

Recall  that  DG{Xo,  uq)  6  £(R  x  Wj’®’,R  x  is  an  isomorphism.  Thus,  the  mapping 

X?G(Ao, tio)  ,  00  is  an  isomorphism  of  (0}  x  Wo'*  into  its  image.  From  (2.7).  we  have 
l{0}x  Wj  ’ 

(4.6)  ^l?G(Ao,  'I)  =  (^(lo).  l?„/’(Ao,  uo)t('). 

By  (4.5)  and  (4.6),  we  conclude  that  D,t,G{Xo,uo)  is  an  isomorphism,  and  this  completes 
the  proof.  □ 

For  u'(A)  defined  by  (4.4),  we  define  the  ‘jump’  /(A)  of  w{X)  at  x  =  xo  by 
J{X)  :=  lim  a(X,x,w{Xy(x))  -  Um  a(A,x,  u;(A)'(x)). 
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From  (4.4),  we  have 


(4.7)  <  f’(A,u;(A)),t;  >=  j^[a{X,z,w{Xy)v'  +  f{X,x,w{X))v]dx  =  v(xo)JCX), 

for  any  v  €  Let  u  :=  ■w{X|^).  Then,  u  is  the  solution  of  Problem  4.2,  and  we  have 

(4.8)  <  f  (Afc,  ti).  Ip  >=  X,  u')v'  +  /(A,  X,  u)ip]dx  =  J(X^), 

for  any  tp  €  with  ^(iq)  =  1. 

O 

We  take  0  €  5/,  with  tp(xo)  =  1  and  fix  it.  Since  (A;,,  u^)  6  M/i,  we  have 


j^[a{Xh,x,  Ufc(^))^'  +  fi^h,  I,  Uhix))^]dx  =  0. 


(4.9) 

From  (4.8)  and  (4.9)  we  get 

(4.10)  J(AO  =  ^  /  [a.(x)(ti'  -  u'J’+V'  +  y5.(x)(u  -  ua)-" V]dx, 

A  J  J 


where 


(4.11) 


ao(i) 

Qfl(x) 

a2(x) 

Mx) 


—  ny(A^,  X, u^(x)),  0q(x)  /y( A^,  X,  ua(x)), 

“  X,u'f^{x)),  ^l(x)  :=  ^/yy(A/i,  X,  U;i(x)), 

=  iayyy(AA,J:,<(2:)  +  ei(w'(a:)-«A(i))),  0  <  ei  <  1, 

=  lfyyyi^k<X,Uk{x)  +  €2mx)  -  Unix))),  0  <  £2  <  1. 


From  the  proof  of  Lemma  4.3  and  lim  ||u  —  UhUj^i.oo  =  0,  there  exists  a  unique  ;  G  Wg 


1.30 


A— 0 


such  that  x(xo)  =  0  and 

(4.12)  y^[ao(x)x'i»'  +  0Q{x)zv]dx  =  j jaoix)ip' v'  +  0Q{x)'>pv]dx ,  'iv  G  I'Fg'^fxo]. 
Note  that  u(xo)  =  ua(2Jo)  =  7-  Hence,  u  —  m  =  0  at  x  =  xq.  Thus,  we  obtain 

(4.13)  j^[ao{x)z'{u  -  UhY  +  0q{x)z{u  -  Uh)]dx 

=  j^[ao{x)tp’{u  -  uh)'  +  0o{x)‘tp(u  -  Ufc)]cix. 

0 

On  the  other  hand,  for  any  u;  G  5^  with  w(xo)  =  0,  we  have 


j^[a{Xh, X, u'(x))w'  +  f{Xh,x,  u(x))ii;]dx  =  0. 
Jja{Xh,x,u'^{x))w'  +  f{XH,x,UH{x))w]dx  =  0. 

Therefore,  we  obtain 

(4.14)  ^  f  [a,(x)iu'  -  +  0,{x)(u  -  Uhy'^^w]dx  =  0. 
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From  (4.10),  (4.13)  and  (4.14),  it  follows  that 

(4.15)  |^(A/,)|  <  ^j[ao{x){z'  -  w'){u'  -  u'^)  +  0o(r){z  -  w)(u  -  Uh)]dx 

!  2  r 

+  \Y.  [at(i)(«'  -  -  ^')  +  -  ^)]di  ■ 

Q 

Now,  let  Zf,  €  Sh  be  the  finite  element  solution  of  z  and  plug  it  into  (4.15).  Then,  we  obtain 
the  following  higher  order  error  estimate  of  lJ(Ah)j; 

(4.16)  |J(A;.)|  <  (|  <  D^FiXn,  uh)(u  -  u^hz-  z^)  >  | 

+  ^1  <  Zh  >  |)(1  +  0(1)). 

Let  us  now  see  the  relationship  between  |A  —  A^]  and  |^(A;,)|.  By  Lemma  4.3,  the  solution 
u;(A)  of  (4.4)  is  differentiable  with  respect  to  A.  Differentiating  (4.7)  with  respect  to  A,  we  show 
that  the  function  A  >->  J(X)  is  and  satisfies 

(4.17)  n(io)^'(A)  =<  DxFi'X,  wi'X))  +  D^F(X,  w{X)){dxw(X)),  v>,  e 

where  5au^(A)  €  Wq’°°  is  the  derivative  of  ui(A)  with  respect  to  A. 

Recall  that  (A,  u)  €  A4o,  DxF{X,u)  ^  0,  u  =  w{X),  and  J(A)  =  0.  We  claim  that  J'(X)  ^  0. 

UJ'(X)  =s  0,  we  would  have  Z?a/’(A, u)+Z?„/’(A,  u)(9ju)  =  0,  and  (5iu)(io)  =  0.  Since  DG(X.  u) 

is  an  isomorphism,  we  have  a  contradiction  that  (0,0)  =  (l,9iu).  Since  lim  |A  -  A/,|  =  0.  we 

/»— 0 

conclude  that  J'{Xh)  ^  0  for  sufficiently  small  k  >  0. 

Thus,  we  have  -JiX^^}  =  J(A)  -  J(Aa)  =  (A  -  A/,)J'(A/,  +  ^(A  -  A^)),  0  <  <  1.  and 

(4.18)  |A  -  Ahl  =  I J'(Afc  +  ^(A  -  Afc))|-'  |7(A,)|, 
for  sufficiently  small  /i  >  0. 

We  would  like  to  replace  the  term  17'(A^  +^(A  -  A/,))|~*  by  some  computable  one.  We  first 
note  that 

^'(Aa+^(A-AO)  =  r{X^)  +  ^iX-X^)J"[X^+MX- X^)) 

=  y'(Afc)(i  +  o(i)),  (o<M<i). 

Hence,  we  just  have  to  approximate  J'{X^^)  instead  of  J'(A/,  +  ^(A  -  A/,)). 

O 

Again,  take  rp  ^  with  tp{xo)  =  1  and  fix  it.  From  (4.17),  we  have 

(4.19)  -/'(A/,)  =<  DxF{Xh,u)  +  £)„f(Ah,  u)(aAu).  ip  >  ■ 
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With  (4.19)  in  our  mind,  we  define  the  ‘approximate  jump’  at  i  =  xq  by 

(4.20)  :=<  DxF{Xh,Uk)+  D„F{Xh,UH){dxUh),i)  > 

0 

where  dxUh  €  5/,[io]  is  the  finite  element  solution  of  the  following  equation: 

O 

(4.21)  <  DuF{XH,^h)(dxUk),VH  >=  -  <  DxF(Xh,Uk),VH  >,  Vt;/,  €  5h[io]. 

Now,  we  would  like  to  estimate  |^'(Afc)  — y)'(A*)|.  First,  we  note  that  ||5au  — =  o(l)i 
because  dxu  €  Wg’**  satisfies  (9iu)(io)  =  0,  and 

<  DuF{Xi^,u)(dxu),v  >=  -  <  DxF{Xh,  u),  v  >,  'iv  €  /fo[xo]. 

Subtracting  (4.20)  from  (4.19),  we  see 

J'M-Jkih)  =  <D^F{XH,u^){dxu-dxUH),rP>  +  <Dl,F(Xh,Uk){u-Uh),rp> 

+  <  Dl^F{Xh,Uh){dxu,u  -  Uh),rp  >  +  higher  order  terms. 

Therefore,  we  obtain 

(4.22)  \J\Xh)-JUX,)\  =  o{1). 

Finally,  we  replace  the  term  u  -  in  (4.16)  by  u  -  u^.  Since  u  =  u;(A),  u  =  u;(A/,),  and 
the  mapping  (A  -  A  +  c)  9  A  •-»  u/(A)  6  Wg’"  is  class,  we  have  ||u  -  <  Co|A  -  A;,! 

with  a  constant  Co  independent  of  h.  Therefore,  we  immediately  obtain  that 

(4.23)  ||ti  -  <  ||u  -  u/,11^1  +  C||u  -  =  ||u  -  +  CCojA  -  X^]. 

Combining  (4.16),  (4.18),  (4.22),  and  (4.23),  we  obtain  the  first  elaborate  error  estimate  of 

|A-Afc|. 

Theorem  4.4.  Suppose  that  Assumption  2.11  holds  for  d  >  3.  Let  (X,u)  €  M  be  such  that 
DxF{X,u)  ^  0.  By  Theorem  2.13  we  can  take  a  nodal  point  xo  €  J  such  that  DG(X,  u)  defined 
by  (2.7)  is  an  isomorphism.  Let  (A^.uj^)  €  Mk  be  the  finite  element  solution  corresponding  to 
(A,  u)  with  u(xo)  =  Uf,(xo).  Then,  we  have  the  following  estimate  o/IA  -  A;,(; 

(4.24)  |A  -  Afcl  <  (l  <  DuF(X^,uh)(u  -  u>,),z  -  >  | 

+  i|  <  D^„F(X^.Uf,)(u  -  Zk  >  |)(1  +0(1)). 

where  z  €  Bq  and  z/,  6  5^  are,  respectively,  the  exact  and  finite  element  solution  of  (4.12)  for 

0 

appropriate  ^  €  5/,  with  ii/(xo)  =  1,  and  Jl,iXf,)  is  the  "approximate  jump'  defined  by  (4.20)  and 
(4.21).  □ 


19 


Now,  let  us  consider  the  second  elaborate  error  estimate  of  |A  —  Ah|. 
Again,  we  consider  Problem  4.1  and  4.1pE-  From  (4.1)  and  (4.2),  we  have 


(4.25)  0  =  (X  —  Xk)  <  DxF^ ,Vh.  >  +  <  DuF^{u  -  Uk),Vh  > 

+  i(A  -  AO'  <  DI,F\vh  >  +(A  -  AO  <  DtF^'iu  -  u^},vh  > 

<  D'„F’*(u  -  Ufc)',t/fc  >  +  higher  order  terms,  6  S/,, 

£ 

where  DxF^  :=  DxFiXh,Uh),  D^F*'  :=  D^F{Xh,u^),  etc. 

We  introduce  the  following  auxiliary  equation:  Find  rj  G  R  and  z  £  Hq  such  that 


(4.26) 


j^[aQ(x)z'v'  +  0Q{x)zv]dx  =  /?  <  6io>^  >-  €  Hq, 


<  DxF*^,z  >=  1, 


where  is  Dirac’s  delta  at  xq.  and  ao,  /?o  are  defined  by  (4.11).  For  the  existence  of  the 
solution  of  (4.26),  we  show  the  following  lemma. 

Lemma  4.5.  For  safficiently  sma.U  h  >  0,  (4.26)  has  a  unique  solution  (r/,  r  )  €  R  x  /fg . 

Proof.  First,  suppose  that  we  have  Case  1,  that  is,  ZAuF'*  :=  Z?„F(A;,,  Uh)  6  CIJ/q,  /f~*)  is 
an  isomorphism.  Then,  we  have  <  DxF'',(DuF^)~^(Sx„)  >=<  Sj-,,  (Z?„F'*)“^(DaF*)  0,  be¬ 

cause  of  the  way  of  taking  the  nodal  point  xq  £  J  (see  the  proof  of  [TBl, Lemma  8.1]).  Therefore. 
(4.26)  has  the  unique  solution  /]  '•=<  Z?aF'‘,(Z?„F'*)“*((5io  )  >“\  z  :=  7?(Z?uF'*)“^(i5^5 ). 

Next,  suppose  that  we  have  Case  2,  that  is,  KerDuF'*  =span{t/^}  and  DxF^  ^  ImZJuF'*. 
There  exists  {9,<t>)  £  R  x  such  that  9DxF^  +  {D^F'')4)  =  and  ^  G  R  is  determined 
uniquely.  We  check  that  ^  ImZ)„F^,  and  hence  ^  ^  0.  If  (5i,  G  ImDuF^,  there  would  exist 
w  €  ffo  such  that  (5i,  =  D^F^'w.  Thus,  from  (TBl, Lemma  8.1],  we  obtain  an  contradiction 
0  5io>  ^  >=<  DuF^w,  xp  >=<  DuF^rp,  xu  >=  0. 

Therefore,  in  this  case,  (4.26)  has  the  unique  solution  t?  :=  0,  x  :=  •^xi'  because 

<  DxF'',z  >=<  DxF^,z  >  +9-^  <  DuF''<P,z>=  <  6x.,z  >=  1.  □ 


Now,  let  us  set  w  :=  u  -  u/,  in  (4.26).  Since  <  (5i„,  u  -  >=  0,  we  obtain 

(4.27)  <  D^F^z,  u  -  u/,,  >=  j^[Qo(x)z'{u'  -  u),)  -f-  Po(x)z{u  -  Uh)]dx  =  0. 
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Since  is  self-adjoint,  it  follows  from  (4.27)  that 

(4.28)  -  <  DuF^{u  -  Uh),Vh  >=<  DuF''(u  -  Uh),z  -  Vh  >,  e 

From  (4.25),  (4.28),  and  plugging  the  finite  element  solution  of  the  equation  (4.26)  into 
Vh  in  (4.25),  we  obtain 

Xh)  <  DxF'',-2h>  =  -<£>„/’'*(«- Ufc),  2  -  rfc  > 

+  ~(A  —  A;,)^  <  Dxxf^t  >  +  (A  —  A/,)  <  D\^F*^{u  —  Uh),  z^  > 

+  r  <  D\^F^{u  -  Zh>  -¥  higher  order  terms, 

and  we  have  proved  the  following  theorem. 

Theorem  4.6.  Suppose  that  Assumption  2.11  holds  for  d  >  3.  Let  (X.u)  6  M  and 
DxF(X,u)  /  0.  By  Theorem  2.13  we  ca.n  taice  a  nodal  point  xq  6  J  so  that  DG(X,u)  defined 
by  (2.7)  is  a,n  isomorphism.  Let  {X^^,u^^)  €  Mh  be  the  Unite  element  solution  corresponding  to 
(A,u)  with  u(io)  =  Ufc(2o). 

Then,  we  have  the  elaborate  error  estimate 

(4.29)  |A  -  Afcl  <  (|  <  DuF{Xh,u^){u  -  Uh),z  --  2/,)  >  | 

+  ^l  <  Dl^f'(XH,UH){u  -  Uhf,ZH  >  l)(l  +  0(1)), 
where  z  and  Zfi  are  the  exact  and  finite  element  solutions  of  the  equation  (4.26),  respectively.  □ 

Remark  4.7.  By  Theorem  4.4  and  4.6,  a  priori  error  estimates  of  |A  —  A^j  are  obtained. 
Since  we  have  a  posteriori  error  estimates  for  l|u-  and  \\z-  2/,(|^^,  and  all  terms  in  (4.24) 

and  (4.29)  ate  computable,  those  estimates  are  a  posteriori  error  estimates  as  well.  The  detail 
of  practical  computation  is  given  in  Section  11.  □ 

5.  Numerical  Examples. 

In  Section  5  we  present  several  numerical  examples  and  discuss  some  points  for  implementation 
of  the  a  posteriori  error  estimates  presented  in  this  paper.  Our  first  example  is  the  following 
simple  one: 

Example  5.1. 

j  u"(z) -t- Au(z)  =  sin(5rt)  In  J:=  (0,1), 

1  u(0)  =  u(l)  =  0.  □ 


21 


The  exact  solution  of  (5.1)  is  u(i)  =  around  A  =  x^. 

A  —  X" 

O 

Let  Ai5  be  the  uniform  mesh  of  7  =  (0,1)  with  15  elements.  Let  5/,  C  Hq  be  the  finite 

element  space  of  piecewise  quadratic  functions  over  A15.  We  use  PITCON  to  compute  the  finite 

0 

element  solutions  (A^,u/,)  6  R  x  5/,  of  (5.1).  Table  5.1  shows  a  part  of  the  output  of  PITCON. 

As  \h  is  getting  close  to  x^,  u  is  getting  bigger  and  the  ‘slope’  is  getting  steeper.  In  Table  5.1 
the  term  |A  —  A/,|  stands  for  the  error  between  the  exact  and  computed  A.  Hence,  |A  —  A/,1  =  0 
means  that  A  is  the  continuation  parameter  at  that  step.  We  see  that  whUe  A  <  7.13828,  A 
is  the  continuation  parameter.  When  A  becomes  greater  than  that  point,  u/i(io),  ^0  or 

xq  :=  8/1  {h  —  1/15)  is  taken  as  the  continuation  parameter. 


Table  5.1:  The  output  of  PITCON:  Example  5.1. 


A/. 

|A-A/,| 

.ffg -error 

Wg’^-error 

0.0 

0.0 

3.67717D-4 

1. 15846  D-3 

1.32361 

0.0 

4.24670D-4 

1.33773D-3 

3.17042 

0.0 

5.41741D-4 

1.70611  D-3 

5.16503 

0.0 

7.71425D-4 

2.42834D-3 

7.13828 

0.0 

1.32875D-3 

4.17803D-3 

7.83691 

2.02465D-5 

1.80565D-3 

5.64451  D-3 

8.85394 

2.29519D-5 

3.59614D-3 

1.12780D-2 

9.38086 

2.43528D-5 

7.44952D-3 

2.34080D-2 

9.50916 

2.46938D-5 

1.00926D-2 

3.17268D-2 

9.68198 

2.51532D-5 

1.93655D-2 

6.08957D-2 

9.75356 

2.53434D-5 

3.12921D-2 

9.83730  D-2 

In  Table  5.1  ‘fi’g-error’ stands  for  ||«— if  A  is  the  continuation  parameter,  and  |A- A/,|  + 
||u  -  if  A  is  not  the  continuation  parameter.  Also,  “’-error’  stands  for  jju  - 

if  A  is  the  continuation  parameter,  and  |A  ~  A^j  -H  ||u  -  u/,||.„i.oo  if  A  is  not  the  continuation 
parameter. 

Now,  we  discuss  how  we  estimate  those  errors.  If  A  is  the  continuation  parameter,  we  surely 
can  use  a  usual  method  to  estimate  l|u  -  We  however  use  the  method  presented  in  this 

paper  to  estimate  ||u  —  u/i|l//i  here  to  show  how  we  implement  our  a  posteriori  error  estimates. 

O 

Let  u/,  €  S/,  be  the  piecewise  quadratic  finite  element  solution  of  (5.1)  over  A15.  We  then 
introduce  another  finite  element  space  Sf,  €  /fg  of  piecewise  polynomials  of  degree  4  over  A15. 
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We  compute  the  finite  element  solution  of  the  following  equation: 

(5.2)  I  i-u'kv'h  +  XuhVh)dx  =  +  sia(Trx)vk)dx,  Vii/,  6  Sh- 

Then,  by  Theorem  3.1,  3.2,  and  Remark  3.8,  and  ||u^((^j,oo  indicate  ||u  -  and 

||ti  -  «^||^^«,  respectively. 

Note  that,  as  stated  in  Remark  3.8,  we  do  not  need  to  solve  the  full  system  of  (5.2).  Instead 
of  that  we  just  solve  tiny  equations  defined  on  each  element  because  we  can  assume  that  the 
values  of  the  finite  element  solutions  of  (5.2)  at  the  nodal  points  are  zero  (see  [BR]). 

As  mentioned  before,  if  A  is  not  the  continuation  parameter,  xq  =  Th  ot  xq  =  8h  (h  =  1/15) 
was  taken  by  PITCON  so  that  Uh(zo)  is  the  continuation  parameter  at  that  step.  In  that  case, 
we  use  the  methods  presented  by  Theorem  3.4,  3.5,  and  3.7.  We  first  approximate  the  exact 
solution  {&,Ui)  ^  R  X  Hq  oi  the  equation 

(5.3)  +  XhU\v)dx  +  9  j  Ukvdx  ~  J {u'|^v’  -  +  sin(Trx)v)dx,  'dv  £  Hq, 

with  £^i(io)  =  0  (see  (3.8)).  We  also  approximate  the  exact  solution  U2  6  Hq  of  the  equation 


(5.4)  +  hU2v)dx  =  jy^v'  -  XhUhv  +  sin(irx)t))(ii, 

with  U2{xq)  —  0  for  any  v  €  Hq  with  v(zo)  =  0  (see  (3.22)). 

Let  {dhfUih)  ^  R  X  Sh  and  U2h  €  Sh  be  the  finite  element  solutions  of  (5.3)  and  (5.4), 
respectively.  By  Theorem  3.4,  3.5,  3.7,  and  Remark  3.8,  we  have  the  following  estimates: 

J  |A  -  Afci  +  ||«  -  <  (|0;.|  +  l|£7i,||^.  )(1  +  0(1)), 

I  |A  -  Afcl  +  ||«  -  <  (l^fcl  +  ||C^iaIU..=o)(1  +  0(1)), 


f  |A-Afc|  +  l|u-Uftl|-  <l|f72,||^.  (1  +  0(1)), 

(5.6)  I 

I  |A  —  A^l  +  ||u  —  <  )|172a||^^oo(1  +  o(l)). 

Again,  to  solve  (5.3)  and  (5.4),  we  take  the  element-by-element  approach,  that  is.  instead  of 
computing  full  systems,  we  solve  tiny  equations  on  each  element.  Note  that  in  the  computation 
of  (5.3),  the  value  9h  is  determined  at  first  by  solving  an  equation  defined  by  (5.3)  on  the  two 
elements  which  contain  iq  as  a  nodal  points.  Then  the  rest  of  Uih  is  computed  using  the  obtained 

In  Table  5.2,  we  summarize  the  results  of  computation.  In  Table  5.2,  ‘A^',  ‘^g-error'  and 
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Table  5.2:  The  estimated  errors  of  |A  —  +  l|u  -  Uh\\:  Example  5.1. 


a 

/fg-error 

^g-est^^^ 

^g‘-est(2) 

WQ'°°-error 

<'"°-est<2) 

0.0 

3.67717D-4 

3.67729D-4 

1.15846D-3 

1.15905D-3 

1.32361 

4.24670D-4 

4.24684D-4 

1. 33773  D-3 

1.33856D-3 

3.17042 

5.41741D-4 

5.41758D-4 

1.70611D-3 

1.70757D-3 

5.16503 

7.71425D-4 

7.71448D-4 

2. 42834  D-3 

2.43154D-3 

7.13828 

1.32875D-3 

1.32877D-3 

4.17803D-3 

4.18819D-3 

7.83691 

1.80565D-3 

1.80641D-3 

1.78547D-3 

5.64451D-3 

5.64856D-3 

5.62764D-3 

8.85394 

3.59614D-3 

3.59692D-3 

3.57331D-3 

1.12780D-2 

1.12864D-2 

1.12628D-2 

9.38086 

7.44952D-3 

7.45042D-3 

7.42542D-3 

2.34080  D-2 

2.34293D-2 

2.34043D-2 

9.50916 

1.00926D-2 

1.00936D-2 

1. 00683  D-2 

3. 17268  D-2 

3.17597D-2 

3.17,343D-2 

9.68198 

1.93655D-2 

1.93667D-2 

1. 93409  D-2 

6.08957D-2 

6.09867D-2 

6.09608D-2 

9.75356 

3.12921D-2 

3.12935D-2 

3.12675D-2 

9. 83730  D-2 

9.85787D-2 

9.85525D-2 

Wo’‘*-error’  are  exactly  same  to  those  in  Table  5.1.  For  A^  less  than  7.5,  ‘/fj-est*'*’  and 
est^^^’  stand  for  and  ||wh||y^,^oo,  respectively,  where  Ufi  g  is  the  finite  element  solution 

of  (5.2). 

For  A^  greater  than  7.5,  ‘fl’g'-est^^^’  and  ‘li'J’°®-est^^^’  are  the  estimated  errors  given  by 
(5.5).  Also,  ‘/fg-est^^^’  and  ‘lVo’’*-est^^^’  are  the  estimated  errors  given  by  (5.6).  We  see  that 
the  estimated  error  /fg-est^^^  and  W^g’^^-est^^^  match  very  well  to  the  corresponding  exact  errors, 
.ffg-est^^^  and  li'g’'*  °-est^^^  are  slightly  underestimated.  We  however  notice  that  those  estimated 
errors  are  very  close  to  the  exact  errors  [ju  -  and  \\u  -  uh||^i.oo,  respectively. 

Now,  let  us  turn  into  a  posteriori  error  estimates  of  |A  -  A/,|.  Suppose  that  A  is  not  the 
continuation  parameter.  In  that  case,  as  stated  before,  either  xo  ■=  7/i  or  ig  :=  Sh  {h  :=  1/15) 

O 

is  taken  by  PITCON.  Let  ^  €  5^  be  a  piecewise  linear  function  such  that  i4’(^o)  =  1.  and 
tp(z,)  =  0  if  X,  ^  zg,  where  z,  are  nodal  points  of  Ajj. 

Then,  we  consider  the  following  equation  (see  (4.12)):  z(zg)  =  0  and 

(5T)  J^{—z'v' +  X^^^v)dx  =  +  Xh^v)dx,  'iv  €  Hq  with  u(zg)  =  0. 

O 

Let  2^  6  be  the  finite  element  solution  of  (5.7)  on  Ajs.  Of  course,  we  can  estimate  the  error 
\\z  -  by  a  usual  method  (see  [BR]). 

O 

Next,  let  Wh  6  Sh  be  the  finite  element  solution  of  the  following  equation  (see  (4.21)): 
^h{xo)  =  0  and 


+  >‘hWhVh)dx  =  -  UfcUfedz, 


O 

Vuh  g  Sh  with  Vh(xo)  =  0. 


24 


We  then  compute  the  ‘approximate  jump’  by  (see  (4.20)) 

J'hi^h)  ■=  +  XkWhip  +  uhi))dx. 

By  Theorem  4.4,  the  error  |A  -  A/,|  is  estimated  by 


(5.8)  |A  -  Afcl  <  |4(A,)|-1||2  -  -  u^ll^.d  +  o(l)). 

On  the  second  method  of  the  a  posteriori  error  estimates  of  |A  — A/i|,  we  consider  the  following 
auxiliary  equation  (see  (4.26)): 

(5.9)  J  u^^zdx  =  I,  and  j  {—z'v' +  Xhzv)dx  =  r]v{xo),  'iv  e  Hq. 

O 

Let  (rjfc,  2fc)  €  R  X  5^  be  the  finite  element  solution  of  (5.9).  To  estimate  the  error  \\z  - 
we  consider  the  finite  element  solution  (jj^,  i/,)  6  R  x  5;,  of  the  following  equation:  j  UhZhdx  =  0 
and 

+  Xhhvh)dx  -  fii^vnixo)  =  7?au/,(xo)  +  ~  ^ii2hVh)dx,  'ivh  €  5^. 

Then,  we  have  the  estimate  \\z  -  <  (lijfcl  +  ||xfcllfri  of  1)).  By  Theorem  4.6,  we  have 

the  error  estimate 

(5.10)  |A  -  A;.|  <  ||r  -  ufcll;,.  !|z  -  (1  +  o(l)). 

In  Table  5.3  we  summarize  the  results  of  computation. 


Table  5.3:  The  estimated  errors  of  |A  -  A^|:  Example  5.1. 


Xh 

|A  -  Ad 

|A  -  AhJ-est^^^ 

|A  -  Ahl-est(2) 

7.83691 

2.02465D-5 

2.09092D-5 

2.09092D-5 

8.85394 

2.29519D-5 

2.36144D-5 

2.36145D-5 

9.38086 

2.43528D-5 

2.50152D-5 

2.50153D-5 

9.50916 

2.46938D-5 

2.53562D-5 

2.53563D-5 

9.68198 

2.51532D-5 

2.58155D-5 

2.58156D-5 

9.75356 

2.53434D-5 

2.60058  D-5 

2. 6005  8  D-5 

In  Table  5.3,  A;,  and  |A  —  A/,|  are  same  to  those  in  Table  5.1.  The  term  ‘|A  -  A;,|-est*"'  stands 
for  the  estimated  error  by  the  first  method  (5.8).  Also,  ‘jA  -  A/,|-est(-^’  stands  for  the  estimated 
error  by  the  second  method  (5.10). 


Now,  let  us  consider  the  second  example: 
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Example  5.2. 

f  u"(x)  =  Ae“<^>  /nJ:=(0,  1), 

(5.11)  { 

[  ti(0)  =  u(l)  =  0.  □ 

The  exact  solution  is,  for  A  >  0, 

A  =  2ii^  cos~^  (|)  ,  u(i)  =  In^cos^  (^)  cos“^^/i(r  —  2)))> 
with  0  <  n  <  TT,  and,  for  A  <  0, 

A  =  — 2/i^cosh~^  (^)  ,  u{x)  ~  In^cosh^  (^)  cosh“^^//(x  -  ^)))7 


with  0  <  n  <  oo.  The  exact  solution  has  a  turning  point  around  A  =  —3.5138307...  {{i  = 
2.39935728...). 

O 

Again,  we  use  the  uniform  mesh  Au  and  the  piecewise  quadratic  finite  element  space  Sh¬ 
in  Table  5.4  we  show  a  part  of  output  of  PITCON  for  the  equation  (5.11).  In  Table  5.4,  the 
meaning  of  the  each  column  is  exactly  same  to  that  of  Table  5.1.  Since  the  exact  solution 
manifold  has  the  turning  point,  the  continuation  parameter  was  changed  three  times  in  this 
case. 


Table  5.4:  The  output  of  PITCON:  Example  5.2. 


Ah 

|A-Ahl 

^Q-error 

W^’°°-enot 

8.40542 

0.0 

1.63866D-3 

7.49987D-3 

5.44008 

0.0 

8.22946D-4 

3.57924D-3 

2.49599 

0.0 

2.20383D-4 

8.91786D-4 

0.55874 

0.0 

1.37612D-5 

5.19053D-5 

-0.31518 

9.05067D-10 

4.98056D-6 

1.79948D-5 

-2.65152 

1.03592D-6 

6.32418D-4 

1.85565D-3 

-3.20484 

2.62501D-6 

1.27568D-3 

3.55051D-3 

-3.51384 

6.80290D-6 

2.80848D-3 

7.86645D-3 

-3.16924 

1.36428D-5 

5.66465D-3 

1.69629D-2 

-2.79878 

1.72583D-5 

7.69351  D-3 

2.41313D-2 

-2.40918 

0.0 

9.92705D-3 

3.22005D-2 

-2.01181 

0.0 

1. 26025  D-2 

4.28233D-2 

-1.65256 

2.38990D-5 

1.56130  D-2 

5.56118D-2 

-0.81887 

2.38261  D-5 

2.73707  D-2 

1.07902D-1 

-0.35494 

1.92353D-5 

4. 4 1004  D-2 

1.82042D-1 

Now,  let  us  discuss  how  we  estimate  the  errors.  As  before,  if  A  is  the  continuation  parameter, 
we  consider  the  following  equation: 

+  ^e'‘''UhVh)dx  =  -  j^(u'hK  +  Ae“'‘i}h)di,  Vi)^  €  Sh- 
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Then,  by  Theorem  3.1,  3.2,  and  Remark  3.8,  (|iih((//^  and  (|u(,(|^^^oo  indicate  ||u  -  and 

l|w  -  “fcllvv‘  .~,  respectively. 

0 

If  A  is  not  the  continuation  parameter  u/,(xo),  Xq  :=  7h  or  xg  :=  Bh  (h  =  1/15)  is  taken 
by  PITCON  as  the  continuation  parameter  again.  Let  (dh,^ih)  6  R  x  S/,  and  [/2h  €  S/,  be  the 
finite  element  solutions  of  the  following  equations; 

+  Afce“'>f7ifci)h)di  +  +  Afce“'“i}h)di,  'ivh  G  S^, 

with  17i/,(xo)  =  0  (see  (3.8)),  and 

jjiU2hVh  +  U2hh)dx  =  -  +  A>,e"'*t/fc)dx,  'ivh  6  5/,  with  i/h(xo)  =  0, 

with  U2kixo)  =  0  (see  (3.22)).  Then,  by  Theorem  3.4,  3.5,  3.7,  and  Remark  3.8,  we  have  the 
estimates  (5.5)  and  (5.6). 

Table  5.5:  The  estimated  errors  of  |A  —  A;i|  +  ||u  —  u/,||:  Example  5.2. 


A* 

iT,} -error 

iyd’‘*-error 

8.40542 

1.63866D-3 

1.63864D-3 

7.49987D-3 

7.51304  D-3 

5.44008 

8.22946D-4 

8.22950D-4 

3.57924D-3 

3.58542D-3 

2.49599 

2.20383D-4 

2.20387D-4 

8.91786D-4 

8.93190D-4 

0.55874 

1.37612D-5 

1.37615D-5 

5.19053D-5 

5.19762D-5 

-0.31518 

4.98056D-6 

4.98073D-6 

4.97980 D-6 

1.79948D-5 

1. 80164  D-5 

1.80155D-5 

-2.65152 

6.32418D-4 

6.32790D-4 

6.31403D-4 

1.85565D-3 

1.85655D-3 

1.85520D-3 

-3.20484 

1.27568D-3 

1. 27705  D-3 

1.27309 D-3 

3.55051D-3 

3.55423  D-3 

3.55014D-3 

-3.51384 

2.80848D-3 

2.81408  D-3 

2.80173D-3 

7.86645D-3 

7.87953  D-3 

7.86662D-3 

-3.16924 

5.66465D-3 

5.68141D-3 

5.65102D-3 

1.69629D-2 

1.70139D-2 

1.69809D-2 

-2.79878 

7.6935  lD-3 

7.71847D-3 

7.67615D-3 

2.41313D-2 

2.42075D-2 

2.41704D-2 

-2.40918 

9.92705D-3 

9.92680D-3 

3.22005D-2 

3.22601  D-2 

-2.01181 

1.26025D-2 

1.26020D-2 

4.28233D-2 

4.29783D-2 

-1.65256 

1.56130D-2 

1. 56603  D-2 

1.55881D-2 

5.56118D-2 

5.59372D-2 

5.58390D-2 

-0.81887 

2.73707D-2 

2.74236D-2 

2.73431D-2 

1.07902D-1 

1. 08632  D-1 

1.08644D-1 

-0.35494 

4.41004D-2 

4.41342D-2 

4.40703  D-2 

1.82042D-1 

1.85414D-1 

1.85044D-1 

In  Table  5.5,  we  summarize  the  results  of  computation.  The  meaning  of  the  each  column  in 
Table  5.5  is  same  to  that  of  Table  5.2. 

O 

To  estimates  of  |A  —  A^l,  we  consider  the  finite  element  solution  of  the  following 

equation  (see  (4.12)):  Zh(xo)  =  0  and 

+  ^h^'‘''ZhVk)dx  =  +  Afce“*t^i;A)(ix,  €  S/,  with  Wh(xo)  =  0. 
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where  ^  €  S;,  is  defined  as  before. 

O 

Next,  let  Wh  €  S/,  be  the  finite  element  solution  of  the  following  equation  (see  (4.21)): 
=  0  and 


J^(Whv'f^  +  Xhe'*''WhVH)dx  =  -  J^e'^'^v^dx,  'iv^eSh  with  Vh(xo)  =  0. 

We  then  compute  the  ‘approximate  jump’  J'|^{Xh)  by  (see  (4.20)) 

J'hi^h)  ••=  +  Xke'^''Wh^  +  e“''i/')dx. 

By  Theorem  4.4,  the  error  |A  —  A;,]  is  estimated  by  (5.8)  as  before. 

For  the  second  method  of  a  posteriori  error  estimates  of  |A  —  A/,|,  we  consider  the  finite 

O 

element  solution  (r?^,^/,)  €  R  x  5/,  of  the  following  equation  (see  (4.26)): 

UhZhdx  =  1,  and  +  Xhe'"'ZhVk)dx  =  Vh^hi^o),  €  S^. 

Let  {rjh,  I;,)  6  R  X  5^  be  the  finite  solution  of  the  following:  J  e^^z^dx  =  0  and 

Jji^h^h  +  Afce“'‘l/,i/fc)(ii  -  ^fcUh(io)  =  Vhhixo)  -  +  A/,e“'‘2;,z;fc)di,  Vi;,  e  S^. 

Then,  we  have  the  estimate  \\z  -  zhWhi  <  (|ji/,|  +  ||4||^p(l  +  o(l)). 

By  Theorem  4.6,  we  obtain  the  estimate  (5.10). 


Table  5.6:  The  estimated  errors  of  |A  -  A;,|:  Example  5.2. 


A/. 

|A-A;.| 

|A  -  Afc|-est<*) 

|A  -  Ah|-est^^^ 

-0.31518 

9.05067D-10 

9.07245D-10 

9.07245D-10 

-2.65152 

1.03592D-6 

1. 03587  D-6 

1.03587D-6 

-3.20484 

2.62501D-6 

2.62472D-6 

2.62472D-6 

-3.51384 

6.80290D-6 

6.80112D-6 

6.801 14D-6 

-3.16924 

1.36428D-5 

1. 36350  D-5 

1.36351D-5 

-2.79878 

1.72583D-5 

1. 72446  D-5 

1.72447D-5 

-1.65256 

2.38990D-5 

2.38564D-5 

2.38567D-5 

-0.81887 

2.38261D-5 

2.37475D-5 

2.37480D-5 

-0.35494 

1.92353D-5 

1.91313D-5 

1.91315D-5 

In  Table  5.6  we  summarize  the  results  of  computation.  In  Table  5.6,  the  meaning  of  the 
each  column  is  same  to  that  of  Table  5.3. 


The  last  example  is  the  following  one  which  is  strongly  nonlinear  (see  [R,pl7]); 
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Example  5.3. 


(5.12) 


+  B{\,u)  =  0, 

u(0)  =  u(l)  =  0, 


in  J  :=  (0, 1), 


where  the  functions  A{y),  B{X,y)  are  defined  as  follows: 


A(y)  :=  arctaii(y/2),  B(X,y)  :=  a(A,y)(l  +^(-6(A,y))  +  b(X,y)T]{a{X,y)), 
a(A,y)  :=  AsiB(y)  +  i/cos(y),  b{X,y)  :=  Acos(y)  -  //sin(y), 

=  {  ! +’27^'*’  t'Mo,  O 

We  set  1/  :=  1.0  in  the  computation. 

The  exact  solution  of  (5.12)  is  not  known.  According  to  the  output  of  PITCON,  as  is  shown 
in  [R],  the  solution  branch  of  (5.12)  is  ‘S-shaped’  and  has  two  turning  points. 

Again,  we  use  the  uniform  mesh  A15  and  the  piecewise  quadratic  finite  element  space. 

We  explain  how  the  errors  are  estimated.  As  before,  if  A  is  the  continuation  parameter,  we 
consider  the  following  equation; 

J^{-Ay(uh)<ikVh  +  ByiX,Uh)uhVh)dx  =  -  j^{-A(u'^)v[  +  B(X,UH)vh)dx,  'ivh  €  S^. 

By  Theorem  3.1,  3.2,  and  Remark  3.8,  ||«fc||jyj  and  ||tih||(yj,oo  indicate  ||u  -  and  ||u  - 

respectively. 

If  A  is  not  the  continuation  parameter  «/i(io).  :=  (h  =  1/15)  is  taken  by  PITCON  as 

the  continuation  parameter.  Let  (5^,  Uih)  €  R  x  5^  and  U2h  €  5/,  be  the  finite  element  solutions 
of  the  following  equations: 

j{-Ay{u'^)U[^v't,  +  By{Xh,UH)UihVh)dx  +  9^  J^BxiX^,Uh)vhdx 
=  -  j{-A{u'^)v'^  +  B(Xh,  Uh)vh)dx,  Vvh  £  Sh, 

with  Uih(xo)  =  0  (see  (3.8)),  and 

j(-Ay{u'^)U2^v'f,  +  By{Xh,UH)U2hVh)dx 

=  -  +  B{X^,UH)vh)dx,  'dvh  £  Sh  with  h(xo)  =  0. 

with  U2h{xo)  =  0  (see  (3.22)).  Then,  by  Theorem  3.4,  3.5,  3.7,  and  Remark  3.8,  we  have  the 
estimates  (5.5)  and  (5.6). 
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To  estimates  of  |A  —  A^],  we  consider  the  finite  element  solution  of  the  following 

equation  (see  (4.12)):  Zhixo)  =  0  and 

j(-Ay{u',,)z'f^v'^  +  By{Xh,Uk)zkVh)dx  =  j +  By{\h,Uh)il>Vh)dx,  * 

O  O 

for  all  €  5/,  with  Uh(xo)  =  0>  where  xp  e  S|^  is  defined  as  before. 

9 

Next,  let  Wh  6  Sk  be  the  finite  element  solution  of  the  following  equation  (see  (4.21)): 
u'fc(xo)  =  0  and 

j^{-Ay{u'f,)w'^v'f,  +  ByiXh,Uh)whVh)dx  =  -  j^Bx(Xh,Uh)vhdx,  >^VheSh  with  i;h(xo)  =  0. 
We  then  compute  the  ‘approximate  jump’  Jj^iX^)  by  (see  (4.20)) 

■=  +  By{Xh,  Uh)whxp  +  Bx{Xh,  Uh)Tij)dx. 

By  Theorem  4.4,  the  error  |A  —  A^j  is  estimated  by 


(5.13)  |A  -  A,|  <  |j;(A,)r>(Ci||z  -  z,\\hx\\u  -  u,|l^.  +  CjHu  -  )(1  +  o(  1 )). 

where  Cj  :=  ||>ly«)|Uoo  and  Cj  :=  ^ll>iy,(w;,)(t/''  -  xOlU*- 

For  the  second  method  of  a  posteriori  error  estimates  of  (A  -  A;,|,  we  consider  the  finite 

O 

element  solution  {vhtZk)  €  R  x  5/,  of  the  following  equation  (see  (4.26)): 


J^Bx(X^,ut^)zhdx  =  1,  and 

j^{-Ay{u’^)z'f,v'^  +  By{X^,UH)zhVh)dx  =  v^vnixo),  'iv  6  Sh- 

Let  {fih,Zh)  €  R  X  5/,  be  the  finite  solution  of  the  following:  j  Bx(Xh,  u/,)i/,(ii  =  0  and 

j^{-Ay{u'^)z'^v'^  +  By(Xh,Uh)zhVh)dx  -  fi,,VH{xo) 

=  Vhhiz^o)  -  j^{-Ay(u'f,)zl,v'^  +  By{Xh,Uf^)zhVk)dx,  G  5/,. 

Then,  we  have  the  estimate  ||x  -  <  (|^fc|  +  j|2j,||//j  )(1  +  o(l)). 

By  Theorem  4.6,  we  obtain  the  estimate 


“  A^l  <  (Cl  j|u  —  uixWffi  ||z  —  ZhWii*  ^2||“  ~  )( ^  +  o(  1 )). 

where  Ci  :=  l|>ly(u;^)||/,oo  and  C2  :=  7pyy(u;^)x^||£«.. 
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Table  5.7;  The  estimated  errors  of  jA  —  A/,|  +  l|u  -  u^^\\  and  |A  —  A/,]:  Example  5. .3. 


.fl'o-est^*^ 

Foi-est(2) 

iyo''“-est(2) 

|A  -  A,,|-est'^> 

|A  —  A;,|-est*“* 

0.19521 

1.14416D-4 

4.57789D-4 

0.57696 

1.09716D-4 

3.80008D-4 

0.97248 

1. 33223  D-4 

4. 14384  D-4 

1.14096 

1.56647D-4 

1.56434D-4 

4.84460D-4 

4.84272D-4 

2.13860D-7 

2.13860D-7 

1.79178 

4.69535D-4 

4.68635D-4 

1. 60921 D-3 

1.60849D-3 

9.00308D-7 

9.00307  D-7 

2.46764 

2.51702D-3 

2.50907D-3 

1.05407D-2 

1.05376  D-2 

7.97236D-6 

7.97234  D-6 

2.99606 

1.24576D-2 

1.23968D-2 

5.72413D-2 

5.72317D-2 

6.19293D-5 

6.19290D-5 

3.24514 

4.20951  D-2 

4.18066D-2 

2.55047D-1 

2. 55964  D-1 

3. 17004  D-4 

3. 17003 D-4 

3.08956 

5.04095  D-2 

3.20935D-1 

1.95627 

1.58211D-2 

7.29373D-2 

1.59592 

1.19638  D-2 

1.19317D-2 

5. 39609  D-2 

5.41812D-2 

3.37883D-5 

3.37878D-5 

1.19898 

1.17190  D-2 

1.16916D-2 

5. 79405  D-2 

5.69798D-2 

2.601 12D-5 

2.60120D-5 

0.99403 

3. 20333  D-2 

3.20370D-2 

1.91055  D-1 

1.91055  D-1 

4.70059D-5 

4.70191D-5 

1.05847 

1. 06782  D-1 

1.0684  8  D-1 

7.26805D-1 

7.27113D-1 

1.97460D-4 

1.97328D-4 

1.29395 

7.73202  D-1 

7.89621D-1 

6. 10406  DO 

6. 43073  DO 

1.12810D-3 

1.11374D-3 

In  Table  5.7,  we  summarize  the  results  of  computation.  According  to  PITCON.  turning 
points  are  at  A  =  3.24513871  and  A  =  0.99403398.  The  meaning  of  the  each  column  in  Table  5.7 
is  same  as  before. 
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